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The Cosmic Microwave Background provides our most ancient image of the 
Universe and our best tool for studying its early evolution. Theories of high 
energy physics predict the formation of various types of topological defects in 
the very early universe, including cosmic texture which would generate hot 
and cold spots in the Cosmic Microwave Background. We show through a 
Bayesian statistical analysis that the most prominent, 5° radius cold spot ob- 
served in all-sky images, which is otherwise hard to explain, is compatible 
with having being caused by a texture. From this model, we constrain the 
fundamental symmetry breaking energy scale to be 0o ~ 8 7 x 10^^ GeV. If 
confirmed, this detection of a cosmic defect will probe physics at energies ex- 
ceeding any conceivable terrestrial experiment. 
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The Cosmic Microwave Background (CMB) radiation was emitted from the hot plasma of the 
early universe roughly 14 billion years ago. All-sky, multi-frequency maps of the CMB sky 
made by the Wilkinson Microwave Anisotropy Probe (WMAP) {1,2) reveal Gaussian temper- 
ature anisotropics of the form expected in standard cosmological scenarios, tracing density 
variations of a part in a hundred thousand in the primordial cosmos (i). However, several 
apparent anomalies in the expected Gaussian, isotropic statistical distribution have also been 
found (4-10). One of the most striking is a large cold spot centred on Galactic coordinates 
b = —57°, I = 209°, with a radius of ^ 5° (10-13). It was detected using the Spherical Mexi- 
can Hat Wavelet (SMHW), an optimal tool for enhancing such features, and has a flat frequency 
spectrum, inconsistent with either Galactic foregrounds or the Sunyaev-Zel'dovich effect (13). 
A conservative estimate of the probability of finding such a feature in Gaussian simulations, 
taking the effect of a posteriori selection into account, is only 1.85%. Several radical explana- 
tions, such as huge voids or an anisotropic cosmology, have already been proposed: many have 
been ruled out by other cosmological observations (14-17). 

Here we consider the possibility that the spot was caused by a cosmic texture (18), a type of 
cosmic defect predicting spots in the CMB (19). Cosmic defects are hypothetical remnants of 
symmetry-breaking phase transitions in the early universe, predicted by certain unified theories 
of elementary particle physics. According to these theories, the different species of elementary 
particle are indistinguishable in the hot early universe. As the universe cools, the symmetry 
between them breaks, in a phase transition analogous to the freezing of water. Just as mis- 
alignments in the crystalline structure of ice lead to defects, misalignments in the symmetry 
breaking in unified theories lead to the formation of cosmic defects (20,21). Breaking a discrete 
symmetry produces domain walls and breaking a circle (or ^7(1)) symmetry produces cosmic 
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strings. Textures form when a simple Lie group (like the special unitary goup, SU{2)) is bro- 
ken. They consist of localized, twisted configurations of fields which collapse and unwind on 
progressively larger scales, a scaling process continuing into the late universe. Each unwind- 
ing texture creates a concentration of stress-energy and a time-varying gravitational potential. 
CMB photons passing through such a region receive a red- or blue-shift, creating a cold or hot 
spot (19) with a magnitude set by the symmetry-breaking energy scale 0o- We have used high 
resolution texture simulations and Bayesian statistical analysis to check if the observed spot is 
consistent with a texture. We conclude that it is, and propose further tests. If confirmed, the 
detection will provide a unique window onto physics at extremely high energies. 

Texture unwinding events may be approximated by a spherically symmetric scaling solution 
(27), on a comoving radius r < kt where k is a fraction of unity and r is the conformal time 
when the texture unwinds. Such events lead to hot and cold spots, with a fractional temperature 
distortion 

AT 1 

^(9) = ±e , (1) 
where 9 is the angle from the center, the amplitude is set by e = Stt^G^q and the scale param- 



eter 9c = 2V2k{1 + z)/ {E{z) Jo and E{z) = ^JQm{1 + zf + Qa, with Qm and Q 



the present day matter and dark energy density parameters and z the redshift of the unwinding 
texture. Since the scaling profile is not valid at large comoving radii r, we truncate Equation© 
beyond its half-maximum by matching its value and derivative to a Gaussian. A scale-invariant 
distribution of spots on the sky is predicted (29), with the number of spots of scale 9c or above, 
Nspoti> 9c) = AnuK^ /{39c)- Here u parameterizes the comoving number density n of un- 
winding textures: dn/dr = ut^^. High-resolution simulations of SU{2) textures yield k ^ 0.1 
and v ^ 2 {30). The corresponding fraction of the sky covered by textures with 9c larger than 
r is fs ~ 0.017. 

We perform a Bayesian analysis of the probability ratio p of two hypotheses given the data 
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(see e.g. (24)). The null hypothesis Hq describes the data D as a statistically homogeneous 
and isotropic Gaussian random field (CMB) plus instrument noise. The alternative hypothesis 
Hi describes the data as CMB plus noise and an additional template, T, given by a cosmic 
texture with parameters e and 6c. The evidence is the average of the likelihood L with respect 
to the prior, 11: Ei = Pr(D\Hi) = J Li{&i\Hi)Il(@i)d&i, where the @i are the parameters 
in hypothesis Hi. This formula naturally incorporates an Occam factor favoring the hypothesis 
with fewer parameters. Our template fitting is performed in a circular area of 20° radius centered 
on 6 = —57°, / = 209°, in the three-year foreground cleaned WMAP combined- frequency map 
(WCM) (11) at w 1° resolution (29). 

The posterior probability ratio p = Pr{Hi\B)/Pr{Ho\'D) = EiPr{Hi)/{EoPr{Hi)) can 
be used to decide between the hypotheses. The alternative hypothesis is favored when p > 1 
and rejected otherwise. The a priori probability ratio for the two models, Pr{Hi)/Pr{HQ), is 
usually set to unity, but since we center the template at an a posteriori selected pixel, we set 
Pr{Hi)/Pr{HQ) to the fraction of sky covered by textures. To compute Ei/Eq, we need the 
likelihood and normalized priors. The likelihood function is just L oc exp (— x^/2), where 

= (D — T)^N~^(D — T) and N is the generalized noise matrix including CMB and 
instrument noise (29). 

As a conservative prior on e, we choose < e < 10~^, the latter value being the upper limit 
imposed by the large scale Cosmic Background Explorer (COBE) satellite measurements (23- 
25). The prior for 6c is obtained by normalizing the number of spots distribution, dNgpot / d6c oc 
6c^ to unity between 6min and 6 max- Photon diffusion would smear out textures smaller than a 
degree or so, hence we set 6jnin — 1°- At large scales textures are unlikely because the sky is 
finite: we set 6 max = 15°. 

We find the probability ratio p a; 2.5, favoring the texture plus Gaussian CMB over the 
Gaussian-only model. The data, the best fit template and their difference are presented in 
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Figure (1). The best fit amplitude and scale are e = 7.7 x 10~^ and 9c — 5.1°. Marginal- 
izing the posterior (i.e likelihood times prior, see Figure (2)), we find dc = 4.9°l2i 
€ — 7.313 ^ 10~^ 95% confidence. The value of e inferred in this way from a single 
extreme event is biased by the detection of signals with high noise, i.e., large Gaussian fluctu- 
ations. To check this we generated 500 all-sky, Gaussian CMB simulations (10,29) and added 
one cold texture spot to each, with amplitude e = 4 x 10~^ below the upper limit, 5 x 10~^ 
inferred from the observed CMB anisotropy spectrum (26). We perform the same template fit 
we applied to the data on each cold spot and then select the spots with high posterior probability 
ratios, p > 1. The mean amplitude obtained from these spots is e ft! 7.9 x 10~^, hence there 
is significant overestimation. Moreover, a more realistic model of textures would predict some 
dispersion in the spot strength, with stronger spots caused by asymmetric, multiple, or moving 
textures. Again, estimating e from the strongest texture spot would lead to a biased value. 

As a complementary test for the a posteriori selection of the template centre, we should also 
study whether prominent Gaussian CMB spots produce such high values of p. Following the 
same procedure for 10,000 Gaussian simulations we select the most prominent spot, finding that 
these spots show typical values of p ft 0.14 < 1, with only ~ 5.8% of the simulations showing 
spots with p > 2.5. Since the kurtosis of the data shows a more significant departure from 
Gaussianity, the percentage may further decrease if spots of all sizes were taken into account. 

In order to analyse further the CMB signal from textures, we generate 10,000 texture plus 
Gaussian CMB and noise simulations, and repeat the analyses performed on Gaussian simu- 
lations with no textures. Considering that we observe one 5.1° texture in about 40% of the 
sky of the WCM - the region unaffected by the Galaxy - we predict around 68 hot and cold 
spots above 9min = 1°- We generate 68 spots per simulation with Nspoti> dc) oc Oc"^ and 
9min — 1°> assigning a random sign and position to each spot on the sky. The amplitude is 
set to e = 4 X 10~^ following the discussion above. We then repeat the previously performed 
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multiscale analysis of skewness and kurtosis (10,11). As there are on average the same number 
of hot and cold spots, the skewness is little affected. However, on the contrary, the kurtosis is 
increased so the anomalously high kurtosis of the data at scales around 5° is actually compatible 
with the Gaussian CMB plus textures interpretation (see Figure 3). 

Fromtheresultsofthetexturesimulations, the predicted number of spots A^5pot(> dc) ~ 11 
for 9c ~ 5.1°, consistent with the single observed spot. Hence we find that the abundance, 
shape, size and amplitude of the spot are consistent with the texture interpretation. The sym- 
metry breaking scale corresponding to the inferred amplitude is 0o ~ 8.7 x lO^^GeV. From 
the relations given below Equation O, the observed texture unwound at z ~ 6, after the reion- 
ization of the intergalactic medium and potentially within reach of very deep galaxy or quasar 
surveys. 

Further observations could test the texture hypothesis. First, if the spot is due to a texture, 
it was caused by a time-dependent gravitational potential. There would be no associated CMB 
polarization. However, if the spot is a rare statistical fluctuation in the primordial density, a 
correlated polarization signal, namely a preference for a radial pattern of CMB polarization 
around it, is expected. On these scales, for adiabatic perturbations with standard recombination, 
almost half the polarization signal is correlated with the temperature anisotropy (27). 

Secondly, there should be many smaller texture spots, with the scale-invariant distribution 
described above. These would be masked or confused by the background Gaussian signal where 
it has maximal power, at ~ 1° — 2°. Nevertheless, each spot would deviate from the expected 
polarization-temperature correlation and a combined all-sky measurement might show a signif- 
icant difference from the standard prediction. 

Finally, a texture at 2 ~ 6 would gravitationally lens objects behind it with a lensing angle 
of order e radians. In particular, it would lens the second-order CMB anisotropics (the Vishniac 
effect) which peak at these angular scales. 
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While certainly radical, the texture hypothesis is the most plausible explanation yet pro- 
posed for the spot. Our analysis shows it to be favored over a purely Gaussian CMB. Alternate 
explanations, such as voids with radius around 200 h~^ Mpc, are far more radical and seem 
inconsistent both with standard cosmology and with galaxy survey observations (28). 
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Figure 1 : A: Azimuthal projection of a 43° x 43° patch of the WCM, centered at (b = -57°, / = 
209°). B: Best fit texture template. C: WCM subtracting the texture template. The temperature 
units shown in the colorbars are jiK and the pixel is 13.7 arcmin. The y-axis points to the 
Galactic north pole. (The template is available on request). 
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Figure 2: Posterior (likelihood times prior) as a function of amplitude e and scale 9c for the 
texture template fit. The prior limit on the amplitude is marked by a dashed white line. 
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Figure 3: Kurtosis of the WCM convolved with the SMHW at 15 scales. The bands represent the 
68% (red), 95% (orange) and 99% (yellow) acceptance intervals given by 10,000 simulations 
of Gaussian CMB (A) and Gaussian CMB plus textures (B). The WCM data (circles) show 
deviation from the expected values compared to Gaussian simulations, but are fully consistent 
with the Gaussian CMB plus textures interpretation at all scales. stands for scale number. 
The 15 scales are: Ri = 13.7, R2 = 25, = 50, R4 = 75, R^ = 100, Re = 150, Rj = 200, 
Rs = 250, Rg = 300, Rio = 400, Ru = 500, R12 = 600, Ris = 750, Ru = 900 and 
i?i5 = 1050 arcmin. 
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Supporting online material 

Supporting Text 
Data and Simulations 

We use the foreground-cleaned WMAP combined map (WCM) of the three year WMAP data 
release (SI). This map is a linear combination of the three maps at frequencies where the CMB 
is the dominant signal, namely 41,61 and 94 GHz. The kpO mask is applied to the data in order 
to exclude foreground contaminated pixels. This mask is the most conservative one of those 
proposed by the WMAP team and excludes around 23% of the sky including a few hundreds of 
known point sources. In the wavelet analyses a set of extended masks (S2,S3) is used in order 
to exclude those pixels convolved with pixels in the mask. The extension is performed around 
the Galactic part of the mask and depends on the wavelet scale R. Here we use an extension of 
2.5R although the results are quite independent of the exact extension (52). Hence for R = 5° 
the mask excludes about 60% of the sky. 

In the template fitting part, we use the WCM in the HEALPix pixelization scheme (S4) with 
resolution parameter Nside = 64. Since the scale of the cold spot we are interested in, is around 
5° (diameter of ~ 10°) this resolution is good enough and reduces the number of pixels used in 
the template fitting. 

Gaussian simulations for each frequency channel are produced, based on the best fit power 
spectrum given by the WMAP team (SI). After convolving with the adecquate beam, uncor- 
rected Gaussian noise realizations are added, taking into account the number of observations 
per pixel. The maps are then combined as for the data. The resolution parameter used in the 
wavelet skewness and kurtosis analyses, is Nside = 256. 
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Template fitting 

We perform a template fitting in a circular area of 20° radius centered at Galactic coordinates 
(b — —57°, I — 209°) in order to check whether there is a texture-like template in the WCM at 
the position of the spot. 

Excluding from the analysis the point sources masked in the three year kpO mask, the total 
number of pixels considered is 1438. Although the angular size of the cold spot is about 10° 
we have to consider at least a 20° radius patch to take into account the whole neighborhood for 
the fit. The Spherical Mexican Hat Wavelet (SMHW) convolves all the pixels in this region and 
they could contribute in an important way to the detected structure. 

To compute the generalised noise matrix of the template fitting, N, which describes the 
correlations between all the pixels, the CMB and the noise contributions have to be worked out. 
The calculation of the latter is straightforward since the number of observations per pixel and 
hence the noise variance is known. As the instrumental noise is uncorrelated the noise matrix is 
diagonal In order to obtain the CMB contribution, we calculate the correlation function for the 
WCM taking into account the pixel and beam effects. As a complementary test we calculate 
the CMB correlation matrix through 70000 Gaussian simulations. Comparing the resulting 
differences of logarithmic evidences obtained from the WCM correlation function with those 
obtained using simulations, the errors are below 1%. 

Choosing different extrapolations of the temperature profile of Equation (2), as an exponen- 
tial or a SMHW, we obtain Aln£^ values between 4.7 and 5.2, which still give p > 1. 

Number of textures 

The expected number of hot and cold spots due to textures, with an angular size larger than 
some 9c (measured in radians), is given by 




(2) 



14 



where the number of texture unwindings per comoving volume, dn, per conformal time, r, is 
(dn/dr) = vt^^, with v a dimensionless constant z/ ^ 2 as measured in numerical simulations. 
Here, Tq is the present conformal time, k a fraction of unity and the factor 47r(ro — r)^ is the 
comoving area of the sphere of currently detected CMB photons at the conformal time r when 
the texture unwinds. If the unwinding event is inside the sphere, the photons "fall in" to an 
outgoing spherical shell of stress-energy and a blue spot is produced. Whereas if the event is 
outside the sphere, the photons "climb out" of the ingoing shell and a red spot is produced. 
The upper limit r < kt is imposed to account for the finite size of the region described by the 
single-texture scaling solution, and the factor two accounts for hot and cold spots. The angular 
scale subtended by the resulting hot or cold spots is given (in the small angle approximation, 
and assuming a spatially flat universe) by ^ = r/(ro — r). If we consider spots larger than 
some size 9c, this imposes the lower limit r > 9c {tq — r) and the r integral is nonzero only for 
conformal times r > tq9c/{k, + 9c). 

We shall be interested in the regime where 9c ^ n, in which case the integral simplifies to 

Nspot ~ -:r^-^- (3) 
6 t>c 

The number of hot and cold spots of angular radius between 9c and 9c + d9c is just the differ- 
ential of the previous Equation, namely 



d9c 3 91' 

which corresponds to Equation (2) in the main text. It follows that the expected fraction of the 

sky covered by spots of angular radius greater than 9c is 

where we approximated the upper limit, where the small angle approximation breaks down, as 
unity. Setting = 2, k = 0.1 md9c = 1° we obtain fs = 0.017. 
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Theories giving rise to texture 

As mentioned in the article, texture is produced whenever a simple Lie group such as SU (2) is 
completely broken. However, there are two important cases to distinguish. The conceptually 
simpler case is that of a global symmetry, in which case the texture evolves as described in the 
article. The order parameter is initially random, with a microscopic correlation length, but it 
becomes aligned on progressively larger scales, with the correlation length growing at a fixed 
fraction of the speed of light. However, if the symmetry is gauged, then the gauge field can 
simply relax to cancel spatial gradients of the order parameter (called a Higgs field in this case), 
so that the texture quickly relaxes to a zero energy density state. The texture winding number 
can still be important - for example, in the standard electroweak theory, it corresponds to the 
baryon number of the universe - but it does not continue to evolve in the late universe. 

It is sometimes argued on general theoretical grounds that the only fundamental symmetries 
of nature should be gauged, not global. For example, it is believed that string theory would never 
predict an exact global symmetry. However, global symmetries are perfectly consistent with 
quantum field theory and many of the simplest theories known do produce texture, even if they 
do not possess any scalar fields. As an example, consider N massless Dirac fermions coupled to 
a Yang-Mills field with a nonabelian gauge group G. This theory possesses a global SU (N) l x 
SU{N)ii chiral symmetry. The full symmetry would be unbroken at high temperature in the 
deconfining phase, but would break to the diagonal subgroup SU (N) as the system cooled 
below the confining transition, leading to the formation of texture. One interesting example 
of a hypothetical global symmetry is family symmetry, for which a natural candidate is SU (3) 
since there are three families of elementary particles. In this case, the SU (3) cannot be gauged 
because then it would be anomalous (S5). Finally, there are known mechanisms in which a 
theory which has only gauge symmetry at a fundamental level can, with exponentially tiny 
corrections, lead to an effective global symmetry (S6). Therefore, while the general theoretical 
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arguments are certainly important, they are not conclusive. Of course, this makes the potential 

observation of a texture even more interesting and significant. 

Collapse and unwinding processes: Nonlinear sigma model 

As the universe expands, the gradient energy in the symmetry-breaking field (or order param- 
eter) is red-shifted away: the initial, random configuration evolves in a scaling manner where 
the field progressively orders itself on a scale set by the Hubble horizon. If the vacuum mani- 
fold possesses a nontrivial topology (for texture, a nontrivial third homotopy group 713), then in 
some regions of space there will be a topological obstruction to field ordering. In such regions, 
the only way for ordering to occur is for the winding configuration to draw itself together and 
collapse down to a microscopic scale, so that the field gradients become strong enough to pull 
the field off the vacuum manifold and over the potential energy barrier. 

The process of field ordering may be described with excellent accuracy by the nonlinear 
sigma model (NLSM) (S7). For the simplest case of SU (2) texture, the space of vacua is a 
three-sphere and the appropriate NLSM is just the 0{N) vector model with 4, with the 
order parameter being a four-component vector whose length is fixed: 0^ = 0o- The equation 
of motion in an expanding universe is 



and this is the equation we solve numerically. Since the textures of interest all unwound in the 
matter-dominated era, long before the dark energy dominated, we use a flat FRW metric with 
scale factor a oc r^. 

The numerical algorithm used to evolve the texture was described in detail in (S8), with one 
minor detail. The "spin flip" procedure described there was, as noted by Borrill et al. (S9), 
unnecessary so it is not used. The code has been parallelized for runs of up to 1024^ on the 
COSMOS supercomputer and is made publicly available on the website 




(6) 
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|http://www.damtp.cam.ac.uk/cosmos/viz/movies/neil.html| In the new simulations textures as 
small as 0. 1 of the horizon are now well-resolved, yielding a considerably higher value of v and 
lower value of k compared to early simulations conducted over a decade and a half ago (SIO). 
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